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RECOGNITION OF BINARY X-RAY SYSTEMS UTILIZING THE 

DOPPLER EFFECT 

B. L. Novak 


Studied herein is the possibility of the recognition of the 
duality of a single class of X-ray systeins utilizing the Doppler 
effect. The procedure is based on the presence of a period which 
coincides with the orbital period at the intensity of the radiation 
in a fixed energy interval of the X-ray consonant of a binary 
system. 

I. We will use 'i'(t,E) to designate the density of the flow 
of radiation from the X-ray source, that is, the average number 
of photons striking a site of unit area, perpendicular to the flow, 
per unit of time in the vicinity of the current time t in a unit 
energy interval in the vicinity of the energy E. One of the 
following forms of the functional dependence of 4' on E is usually 
adopted [1] : 


mE)= .■ 


Let 4i(t,E) and L(t), respectively, be the energy and time 
spectra, that is, the average number of photons recorded by the 
radiation detector: 

-in a unit energy interval per unit of time for the energy 
spectrum, 

-in a unit of time for the time spectrum. 

The energy spectrum is associated with the density of the flow 
by the relationship 


*Numbers in the margin indicate pagination in the foreign text. 
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Zi 


<P(t,E)=£S?(t;E). 


where c is the effectiveness of registration, and S is the area 
of the detector. We will assume that the detector does not rotate 
in absolute space and 


E 


£^=const^ 0 

0 


with 


with 





where and E|^ are some energy values which are characteristic 
for the given detector. It follows from (4) that 4>(t,E) differs 
from zero only with Ej^<E<Ej^. 


It is clear that 

L(t)=Iop(t.£)l!: 


l#(t,EUE. 


E. 



We will use n=n (t^ , t^+At, E{ , E^ ) to designate the calculation 
of the pulses with energies from EJ to E| in the detector within 
the time interval of a magnitude At from t^ to t^^+At, where 
t^=t^+ (i-1) At , and t^ is the beginning of the observation process. 
Then, the process of registration of the radiation may be described 
in the following manner: the sequence of measurements {h^} is gener- 
ated by the Poisson process with variable intensity 



Here, G(E',E) is the function of response of the detector, 
that is, the conditional density of distribution of the pulse 
energy, occurring during registration of a photon of a certain 
energy by the detector. 
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2. If the X-ray source is part of a binary system, then its 
flow density 'f(t,E), and, consequently, its energy spectrum (t>(t,E), 
changes with time for two reasons: 

1) the orbital motion of the source, which evokes the Doppler 
effect, 

2) the variability of the radiation process itself. 

We will first examine the form of the energy spectrum (|)^(t,E), 
the variability of vdiich is associated with the Doppler effect. ^ 

Let <I'q(E) be the "unperturbed" spectrum, that is, the spectrum 
which would have a source if it were not part of a binary system. 

't is not difficult to show that 



Here, v(t) is the projection onto the line of sight of the velocity 
of the X-ray component, relative to the center of mass of the 
binary system, and c is the magnitude of the speed of light. 

We would note that 

-jtfflJU (») 

o 

and does not depend on t. 

We have indicated three forms of the functional dependence of 
density of the radiation flow on the energy. Three forms of 
spectra 4 j (t,E) correspond to these forms of the dependence. 


the 

the 
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If the density of the radiation flow is described by expression 
(1) , then 


and with v(t)<<c 

If the density of the flow is described by (2), then 


and with v (t) < <c 






E))E^ 


r 



If the density of the flow is described by (3), then 

tp, tt, Q = a* E 
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emd with v(t) <<c 


We will now switch to the energy spectrxun which describes 
the variability of the radiation process. We will designate this 
spectrum using 4>^(t,E). According to [2-5], there exists a class 
of X-ray sources in vdiich the variability of the radiation process 
has a random nature. It is precisely this class with which we 
will be occupied subsequently. Judging from those same studies 
[2-5] , the time spectra of the sources of this class are described 
well by a model of the random pulse process {shot noise) , that 
is 




where £(t) are pulses of a given form, and the times t^^ are random 
magnitudes, distributed according to Poisson's law with some 
parameter u . 

If the time spectrum is described by a model of the random 
pulse process, then, keeping relationship (5) in mind, it seems 
natural to utilize a similar model for the description of the 
energy spectrum, as well. We will assume that 


cp,(t,E)=2;:if(t-t,,E). (13) 


where 0(t-t.,E), similar to (12), are pulses of a given form, and 
they have the very same meaning as in (12) . Following study [5] /J_ 

in the assumption of the nondependence of the times of appearance 
of the pulse 4»(*) on the energy, we will assume that ♦(•) can be 
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represented in the form 


<f(t-tc,E)='fi(t-tt)'fJE). ii4) 


We will now construct a general tiodel of the energy spectrum 
^ (t,E) , We will assume that it consists of two additive com- 
£K>nents, associated with the above-indicated reasons for the 
variability. Thus, in the case when the source is part of a 
binary system, 

and in the opposite case 

where ^ is the sign of the mathematical expectancy. 

For the evaluations which we will carry out subsequently, it 
is desirable to know the connection between the characteristics of 
the time and energy spectra. In this connection, we will give our 
model in detail. We woul note that, in virtue of (5), (12), (14) 
and (15) , and (16) , 

and it is natural to assume that £(t) = ^i(t), that is, 

tlth Jj(t-ti,E)JE. 
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Ve will asavsne that the distribution of the pulse i ( • ) ^ according 
to the energies, takes place in accordance with the energy spectrun. 
0 (E) , that is, 

J4UEME . 



or 


a 


It is evident that the model we constructed generalizes the 
model in [2] in a natural manner. In principle, one can compli- 
cate the model by taking into account the dependence of the dis- 
tribution of the pulse £(•) according to the energies on the 
true spectrum $(t,E), and not on (E) . 




3. We will examine the probability properties of the intro- 
duced magnitudes. We would note that, since relationships (15) 
and (16) differ only in their determined additive terms, the 
probability characteristics of the spread for these cases coincide. 
Thus, if <j)(t,E) and <>(t,U) are functions determined by expressions 
(15) or (16) , with various values of the energy E and U, then, 
making use of the results in [2J , we obtain an expression for the 
mutual spectral density S(v,E,U) of the random processes i}i(t,E), 
Mt,U) 


s(v,E,u)=^ifwr 


ai(E)W 

( !We)Je)* 
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Here, P(v) is the transform of the Fourier function i(t). The 
mutual covariation function of the processes ^(t,E) and ^(t,U) is 
stationary, and is determined by the expression 



The average values of the magnitudes determined fron relation* 

ships (15) and (16) will be distinguished using the indexes g and 

o, respectively. Thus, the mathematical expectancy of the random 

process 4i(t,E) is equal to ^ (t,E) . 

8*0 

From the obtained results and expression (6), we conclude 
that the intensity A (t) is a random process with a mathematical 
expectancy 


K Eh 





and a stationary covariation function 



We will now exeunine the probability properties of the sequence 
of measurements {nj. It is evident that 



We will obtain expressions for the covariation of the magnitudes 
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UKKJtNAL PAGE Ib 
OF FHX>R QlMU'rY 


and Initially, let j>0, that is, the time intervals in 
the course of which the calculation is recorded do not intersect. 
It is not difficult to see that, for this ease. 






Here, cov is the sign of covariation. 

We will make a substitution in the integral of the variables 
Ti*ti; T«t 2 -ti. Then, 

(j-iWt 


'I ti.j 

((.-iUt 

g-iVrt 


( 23 ) 




where 




(t-jat) with T* ^ jat 

with T * 








It is evident from expression (23) that the covariation is 
stationary, that is, it depends only on the magnitude of j. Ve 
will use Cj to designate cov (n^ ,n^^^ ) . 

Now, let j*0. Utilizing the results given in (6], we obtain 
the following express ion for the dispersion of the magnitude nj^: 

i^*et At 

cf =cov(nw,»ij= j M j,c/t)Jt^JfAt-t)8(i)ae. (24) 

■kn o 




spectral density of the seqiMnce of neasur^Mnts {n^} 
is expressed via the covariation and dispersion according to the 
formula 



where c_j*Cj. Specifically, 


t>At 


r5,o(v) = E Cj = I . 


with 



r(v)4r(c). 



4. We will now switch to the verification cf the duality of 
the source. In the current study, for simplicity, we will examine 
only circular orbits. In this case, v(t), which is part of 
relationship (7), is a harmonic function, and Mg(t) for a binary 
source, as follows from (20) and (9) or (10), or (11), may be 
represented in the form 

Uo+aiCOS(0,t+<f)* (28) 

where is the circular frequency of rotation along the orbit, and 
a^,ai, ^are some constants. 

We will examine the sample spectrum of the sequence of measurements, 
that is [7] , 
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Here, M is the total nu^>er of measurements, taking part in the 

processing, » and W(K) is the correlation window. In 

our analysis, we will limit ourselves to the simplest case, when 
W(K)«1 with K>^:1,2, . . . ,M, that is, to the study of the unleveled 
spectrum C (v) [7] . 

Let us assvune that 

gn=|}^Snj , 

and 

C i(>)= ^ ^ '^j - i ^ 

I c^(v)= ^ +2S VHCos2nvK/il,x 

X S(^h^“AS)(An-+K-Ah)) * 

: CjW= ^ CS(|hj-| R)(A*i^-An) +2 EWk<«4,vkA 

Then 

•)R'k'^:aL PA.GE b 
;,p n-v iR O'VALH^ 

.11 

I 


CCv)=c^W+Qfv)^C3Cv). C») 

We will occupy ourselves with the study of Cj (v) . We will assume 
that 

' (so) 


and, for convenience, we will examine the continuous analog Ci (v) , 
that is, 

45 cosa’jnfu. . (3i) 

_ O O - ‘ 

Since relationship (22) is fulfilled, then 



T t + At 


S = T i ^ i . 


/1 2 



For evaluation, we will limit ourselves to the main loop of Ci (u) . 
Then, making use of relationships (22) and (30) , and having 
carried out the operations determined by expressions (31), (32), 
and (33) , we obtain, in the case when the source is binary: 

Max Q(v) ~ At? (34) 


We will switch to the study of C 2 (v). We will evaluate the 
maximally possible magnitude of C 2 (v) . Since C 2 (v) is a random 
function of the argument v, then we will make use of the traditional 
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I 

I 

■1 


probability evaluation 
possible values of this function. Utilizing the results of [7], 
we obtain 

gqw=rtv)- 

amd 

i qcv) + 3 =& + 3\fT) . (35) 

From (35) and (27), the following evaluation follows; 

Q(^)^(d+3)^^)^(o). (=6) 

On the other hand, according to the properties of the sample 
spectrum, C2(v)»0. Thus, 

0 (i+3\Z2r)r(o) 


gQw)*3ykiv)' for the maximally 


and, in that case when the X-ray source is part of a binary 
system with sufficiently large values of T (greater than some T*) , 
the maximum value of the function C(v) exceeds the threshold value, 
determined by the expression (l+S/?) (0) . If the source is 

single, then the values of the sample spectrum are no greater than 
(1+3/7) r (0). 

O 

5. We will evaluate the magnitude of T*, having adopted the 
Gaussian form of the pulse of the time spectrum. Thus, let 

where a and b are some constants. 
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The transform of the Fourier function i{t) is then determined by 
the expression 


^tv) =• afiVFe 




Let the function of response of the detector be ideal, that is, 


g-(e;e)=S(e-e), 


where 6(*) is Dirac's delta- function. Then, the covariation 
fxmction of intensity X (t) , according to (21), is determined by 
the expression 


'(0,(e)(Se] 


Hence, utilizing (26), we obtain C' 

Qo)=J MJrtIt + - . 

^ tj ^ c«o(E)<SEj 


We would note that, according to [2], 


-1! 

felt 


f: Llt)=]iU(t)<St=|Jjxe 


On the other hand. 


i.(t)=6 rn>(V)<IE- i!tUEME. 

Eh 


Therefore , 


0 






For definiteness, let the dependence of ¥ on E be described 
by relationship (1) ; then 





Here, is the maximum value of the projection of the velocity 
of the X-ray component, relative to the center of mass of the 
binary system, on the line of sight. 


Hence, 

fviaxt»=2afTAt 




and T*, according to (34) and (36), is determined from the equality /14 
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We will utilize the obtained formula for the evaluation of T* 
during the registration of the radiation by the detector with the 
parameters of the SNEG-2MP instrument (81 ; 

= E'=20k.v; Ei=280k.r . 

We will take the parameters of the Sud X-1 as the characteristics 
of the X-ray source (5,9] ; 

y =i,S3; Vo=58l'^sec . 

Through the utilization of these data, the value of T* is 
obtained as equal to 4.4 days, which should evidently be recognized 
as acceptable. 
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